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A theory that describes how to load negative charge into a nonlinear, three-dimensional plasma 
Wakefield is presented. In this regime, a laser or an electron beam blows out the plasma electrons 
and creates a nearly spherical ion channel, which is modified by the presence of the beam load. 
Analytical solutions for the fields and the shape of the ion channel are derived. It is shown that very 
high beam-loading efficiency can be achieved, while the energy spread of the bunch is conserved. 
The theoretical results are verified with the Particle-In-Cell code OSIRIS. 



Plasma-based acceleration relies on an underdense 
plasma to transfer the energy from a laser beam or an 
electron beam to a trailing bunch of electrons or positrons 
[TJ 12]. The beam load is accelerated in a wake moving 
with a velocity near the speed of light, c, until the driver's 
energy is exhausted or — in the case of a laser driver — 
until it outruns the plasma wave. Major advances in 
both laser- and beam-driven accelerators have recently 
been achieved in a regime in which the fields of the driver 
are so intense, they push all plasma electrons aside, gen- 
erating a pure ion channel [3]. 

For the laser-driven accelerator [T] , experiments [H |S] 
have inferred and simulations O [7l [8] have demonstrated 
that electrons self-injected into the ion channel can form 
a quasi-monoenergetic beam. Externally injected, low- 
charge bunches have been shown to improve the repro- 
ducibility and the quality of the final electron beam [9]. 
In the beam-driven case stable acceleration of the 
tail of a A2GeV electron beam culminated in the dou- 
bling of the energy of some electrons in less than 1 me- 
ter jlOj. The high-gradient acceleration of a short trail- 
ing electron bunch in the wake of a driving beam, which 
is central to the afterburner concept JT], has also been 
achieved [12]. A theory that describes the wakefield in 
this blowout regime has recently been developed [T31IT4]. 

While there has been tremendous progress both exper- 
imentally and theoretically on understanding how wakes 
are excited in the 3D nonlinear regime, there has been lit- 
tle work on how the trailing beam loads the wake. In the 
linear regime, the issue of beam loading was addressed in 
Ref. [15], where the wakefield generated by the trailing 
bunch was superimposed on that of the driver to yield 
the final accelerating field. Thus, the maximum charge 
that can be loaded was evaluated and the current pro- 
file that makes the wakefield within the bunch flat was 
determined. Ref. [IH] also discussed the effects of trans- 
verse beam loading, emittance, and phase slippage, and 
concluded that beams with spot sizes much smaller than 
those of the wake are required. 



When blowout occurs, the accelerating fleld is identical 
within each transverse slice of the ion channel [T3J [TJ , so 
as opposed to the linear regime, transverse beam loading 
does not affect the energy spread of the beam. Addi- 
tionally, because the focusing force in the ion channel is 
linear, the emittance is conserved. Therefore, the most 
important consideration for reducing the energy spread 
is keeping the accelerating field constant along the prop- 
agation direction. We note that for high-energy physics 
applications, narrow trailing bunches are still needed for 
matched beams and for reducing synchrotron radiation 
losses [TS]. 

In an estimate offered in Ref. [16j , the number of par- 
ticles was found to scale with the normalized volume of 
the bubble (or the square root of the laser power) . The 
same scaling was obtained in Ref. [17] but the coeffi- 
cient, determined by simulations, was more than three 
times larger than the one estimated in Ref. [16]. These 
results are not necessarily contradictory because, in prin- 
ciple, one can choose to accelerate either a small number 
of particles to high energy or a large number of parti- 
cles to low energy. The question of merit is not just how 
many electrons can be loaded, but what kind of electron 
bunch can most efficiently convert the energy available 
in the wake of the driver into kinetic energy uniformly 
distributed to its electrons. 

In Ref. J3], the wakefield in each transverse slice 
was found to be proportional to the product of the local 
radius of the ion channel rf, and the slope drb/d^, where 
(^,rb) are cylindrical coordinates with ^ = ci — z and 
the driver is moving toward positive z. The shape of the 
bubble is represented by the trajectory of the innermost 
particle given by Eq. (11) in Ref. [T3]. This description 
is valid between the points where the particle trajectories 
cross, at the very front and the very back of the bubble. 
To make progress analytically, we take the ultrarelativis- 
tic limit, where the normalized maximum radius of the 
ion channel is ujpRi,/c 3> 1. The equation for the inner- 



2 




FIG. 1: The electron density from a PIC simulation with 
OSIRIS [18] for kpRb = 5 is presented. The beams move to 
the right. The broken black line traces the blowout radius in 
the absence of the load. On the bottom, the red (black) line 
is the lineout of the wakefield Ez{£,,ri, — 0) when the beam 
load is present (absent). 



most particle trajectory reduces to (see Ref. [T3]): 



rb 



(1) 



charge profile of the load is described below. 

If the repelling force is too large and the beam too long, 
the electrons in the sheath will reverse the direction of 
their transverse velocity at some ^r, where ^1^=^^ = 0, 
and consequently £^z(Cr) = 0- This is a very undesir- 
able configuration because it implies that the front of 
the bunch feels a much stronger accelerating force than 
the back. 

We are interested in trajectories for which r\j{^ > 0) 
decreases monotonically. A may then be expressed as a 
function of r?,: A(^) — /(?"(,). Substituting r[,' — 
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where the prime denotes differentiation with respect to 



Eq. ( 1 ) reduces to ^ 
be integrated to yield: 
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Before further analyzing Eq. ([2]), we comment on 
salient features of the unloaded case {l{rh) = 0). Evalu- 
ating the constant in Eq. ^ from the condition Ez{rh = 
Rb) — 0, we obtain: 
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l,Rb>rb>0 (3) 



Eq. ^ can be integrated from the top of the bubble 
'''biC — 0) = Rb to yield the innermost particle trajectory 
for < Tf, < Rb'. 



where we adopt normalized units, with length normal- 
ized to the skin-depth c/cjp, density to the plasma den- 
sity Up, charge to the electron charge e, and fields to 
mcWp/e. The term on the right hand side of Eq. ([T|) can 
describe the charge per unit length of an electron beam 
and/or the ponderomotive force of a laser [13j. Here we 
are interested in the back half of the bubble, where the 
wakefield is accelerating and the quantity 27rA(^), with 
A(^) = /p rubdr, is the charge per unit length of the 
beam load. 

We define ^ = at the location where rj, is maximum, 
i.e., ^1^=0 = 0. In Ref. [T^, it was shown that for 
ujpRb/c ^ 1, the wakefield is ~ therefore, 
Ez{^ = 0) ~ 0. For ^ > 0, the electrons are attracted by 
the ion channel back toward the ^-axis with ^|{>o < 
until ^ = where beam loading starts. For S, > ^s, 
the electrons feel the repelling force from the charge of 
the accelerating beam, in addition to the force from the 
ion channel. The additional repelling force decreases the 
slope of the sheath ^ , thereby lowering the magnitude 
of Ez- This can be seen in the simulation results in Fig. 
[l] where the trajectory of the innermost electron for an 
unloaded wake is drawn on top of the electron density for 
a loaded wake, and the corresponding wakefield for the 
two cases is also plotted. The method for choosing the 
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where F{ip\m), E(ip\m) are the incomplete elliptic inte- 
grals of the first and second kind [l^ . 

To minimize the energy spread on the beam, we 
seek the beam profile that results in E^ifb < rg) — 
h^b^\ri,=rs — constant = —Eg within the bunch. The 
shape of the bubble in this case is described by the 
parabola = - 4Esi^ - 6)- For < ^ < E^ 
is given by Eq. ([s]). Eg is found by requiring that the 



wakefield is continuous at S,s'. Eg 
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< ^ < -I- , where £,s + is the location at which 
the sheath reaches the ^-axis, the profile of A(^) that leads 
to a constant wakefield is trapezoidal with maximum at 



Ej 



and minimum at A(^s 
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and he total charge Qg = 2tt J^^^^'^^"'' X[^)d^ is: 
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Eq. ([6| illustrates the tradeoff between the number of 
particles that can be accelerated and the accelerating gra- 
dient. Because Qs x Eg is the energy absorbed per unit 
length, Eq. ([6]) also indicates that the efficiency from the 
wake to the beam load does not depend on the field Eg, 
as long as the charge profile is chosen appropriately. In 
Fig. [T| the charge of the beam load was chosen using Eq. 
^ for Es — 0.35Rb- For this E^, the location can be 
obtained either from a simulation for an unloaded wake 
or from Eq. (|3|-(|4|, and the charge profile from Eq. ([5|. 

In a linear wake, a wide electron bunch with total 
charge Qi and transverse spot size A, loaded at some 
^0 with the appropriate profile [15] can also lead to a 



flat Wakefield (Ei 



e no 



cos(fcp^o)) within the bunch. 



The total accelerating force is: Qi x Ei — — ^), 

where g;^^ is the energy per unit length of the wake in 

front of the bunch, and -^A that behind it. The effi- 

ciency, (1 — increases for a decreasing accelerating 

gradient and reaches 100% for Ei = 0. This is in stark 
contrast to the blowout regime, where the efficiency rij, is 
constant for any E^. 

To calculate the efficiency ry^, let us assume that the 
bunch is terminated at some where A^f = — < 
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After this point the wakefield is described by Eq. 
(pi) with l{rh < Tf = rb{£,f)) = 0. From the boundary 
condition Ez{rf) = —Eg, we obtain Ez{Q < rf, < rj) ~ 
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In 



Rcf. [16 , it was shown that the energy in the fields of a 
bubble is proportional to R^. Therefore, the energy per 
unit length available to the bunch Sq;^^^! scales as Rf, 
and that left behind it Siq^^ scales as R^. We define the 
beam-loading efliciency as rib =Jf avail " '^lost)/'^avail- 
Substituting the expression for Rb'. 
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where Qg is the charge of a trapezoidal bunch that is de- 
scribed by Eq. (l5| but terminated at -f j instead of 

2 — 

^« + We note that the efficiency approaches 100% 

2 — 

for A^y jg- ^ Qs Qs- Because the mathematical 
formulation involves approximations, there is still some 
energy in the plasma behind the bunch, even with the 
optimal A(^). This is the case in Fig. [l] where a second 
bubble with radius Rb2 ~ Rb/"^ does appear, but because 
QsX Eg (X Rf, the efliciency is still rjb ~ 90%. The wake- 
fleld within the bunch in Fig. [l]is constant, in agreement 
with the theory. We note that a 10% deviation of the 
total charge for a flxed bunch length leads to a wake that 
is no longer flat. 

It is illustrative to compare the amount of charge that 
can be loaded into linear and nonlinear wakes. If we 



assume for the linear wake an effective A ~ (? juJ^, which 
is required for high efficiency and good beam quality [15] , 
we have: 
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where = e^/(mc^) is the classical electron radius. In 
the linear regime, the density perturbation is rii/riQ <C 
1. In the blowout regime, because the total accelerating 
force scales with the fourth power of the blowout radius, 
a radius kpRb ^ 5 leads to a total force ~ 1000 times 
larger than that in the linear regime. Eq. ^ can be 
converted into an engineering formula: 



Qs ^ eEs 
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For a bi-Gaussian beam driver with fcpcr^ ~ 1 and 
kpffr < 1 we have kpRb ~ 2^^(fcp(Tr)^ = 2V7V [Ti] . 

and for a matched laser driver kpRb ~ [14| . where 

flo is the normalized vector potential. Both for a beam 
driver with 3 x 10^° electrons and CTr- <C cr^ = 16.8/im 
in a plasma with rip — lO^^cm"'^, and for a matched 
laser-driver with power P = 200TW^ in a plasma with 
rip — 1.2 X lO^^cm^'^, we have kpRb ~ 4. Choosing 
_eE^ = ~ 2, we obtain Qs ~ 1.9nC for the beam- 
driven case and Qs — O.bbnC for the laser-driven one. 

Another analytically tractable case is that of a bunch 
with a flat-top proflle starting at ^ = ^j: Z(0 < < 
rj) = Iq. The behavior of the flat-top bunch is important 
because it is similar to that of a Gaussian bunch (see 
below). For such bunches, Eq. ([2| becomes: 



E,{0 <rb< rs) 
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There are three distinct cases, all of which can be 
solved analytically. For small charge per unit length 
Iq < i?^/(8r|), the plasma electrons reach the ^- 
axis quickly with some remaining kinetic energy. For 
Iq > i?^/(8r?), there is a minimum radius r,„ = 



4^0 - \/(4'o)2 + - Slorl for which E^{^rn) = and 
the transverse velocity of the innermost particle changes 
sign as described earlier. At the bunch must be ter- 
minated; otherwise, for ^ > ■^m, it will experience a de- 
celerating fleld. Because the plasma electrons do not re- 
turn on the ^-axis, they still have some potential energy. 
Thus, for Iq ^ i?^/(8r?), there is always some energy in 
the plasma behind the bunch. 

For a flat-top bunch, the beam-loading efliciency is 
maximized if Iq = i?j/(8r?). Then the shape of the bub- 
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ble and the wakefield are given by: 
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and the innermost particle will reach the ^-axis at + 
A^s, where A^j =^{Rl- ^Rf-rl). In this case, the 
energy absorption per unit length is identical to that of 
an optimal trapezoidal bunch 27r/oA^j(|£'2|) = Qs x Es- 
The difference in the accelerating force experienced by 
the front and the back of the bunch will tend to increase 
the bunch's energy spread. This can be avoided either 
by injecting the bunch with an initial energy chirp to 
compensate for the effect caused by the field in Eq. (13) 
or by using a monoenergetic trapezoidal bunch. 

If the driver travels with a velocity slower than that of 
the accelerating electrons, these electrons will move with 
respect to the wake. In this context, it is interesting to 
see what happens if a flat-top electron bunch optimized 
for some ^1 is instead placed at ^2 and ^3, both smaller 
than ^1. 




FIG. 2: Wakefield lineouts for (a) a flat-top electron bunch 
and (b) a Gaussian bunch with the same charge at three difi'er- 
ent locations ^1 (red), ^2 (blue), ^3 (green) is plotted from the- 
ory (solid lines (a)) and simulations (symbols (a), (b)). 



In Fig. [2f^a), we compare the lineouts of the wakefield 
Ez{^, Tb — 0) from three 2D cylindrically symmetric sim- 
ulations with the theoretical results for flat-top beams. 
For each simulation, an electron bunch with /q = 0.25i?^ 
and length A^j = 0.27i?b is loaded at one of three loca- 
tions: ^1 = 0.67i?h,6 = 0.53i?b,C3 = 0.31i?6. The open 



red squares correspond to loading at ^1, the solid blue di- 
amonds to ^2, and the open green circles to ^3. The solid 
lines are derived from the theory (for Iq > i?^/(8r?), the 
particle trajectory in the region ^ ^ < ?m can be writ- 
ten in terms of the integral E{ip\m)) and are in excellent 
agreement with the simulations in all three cases. 

We repeated the simulations using Gaussian bunches 
with the same number of particles as in the flat-top cases 



and iVb(z) 



/2ir cr^ 



where = A£,s/{2^/2). 



Each bunch is placed so that its center is at a distance 
V^fTz from ^1,^2, and ^3 for the three simulations. The 
results, shown in Fig. |2]^b), conflrm that the Gaussian 
bunches may be treated using the theory for flat-top 
bunches. In both Fig. [2f^a) and[2jb), we observe that 
the wakefield is relatively flat regardless of the placement 
of the bunch. The initial negative slope is balanced by 
a smaller positive slope for most of the acceleration pro- 
cess. 



The main approximation in the theory, kpR^ ^ 1, is 
expected to hold for kpRb ^ 3. However, the formal- 
ism described here can still be applied if one numerically 
solves Eq. (11) of Ref. [13]. 
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